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Preface 



Because of the immense interest in and success of the first edition, the second edition of Advanced 
Engineering Electromagnetics has maintained all the attractive features of the first edition. This 
edition contains many new features and additions, in particular: 

• A new chapter, Chapter 14, on diffraction by a wedge with impedance surfaces 

• A section on double negative (DNG) metamaterials (Section 5.7) 

• A section on artificial impedance surfaces (AIS, EBG, PBG, HIS, AMC, PMC) (Section 8.8) 

• Additional smaller inserts throughout the book 

• New figures, photos, and tables 

• Additional examples and numerous end-of-chapter problems 

Purchase of this book also provides you with access to a password-protected website that contains 
supplemental multimedia resources. Open the sealed envelope attached to the book, go to the 
URL below and, when prompted, enter the unique code printed on the registration card: 

[http://placeholder.for.actual.url.tk.com] 

Multimedia material include: 

• PowerPoint view graphs in multicolor, over 4,200, of lecture notes for each of the fifteen 
chapters 

• Forty-eight MATLAB® computer programs (most of them new; the four Fortran programs 
from the first edition were translated to MATLAB®) 

Given the space limitations, the added material supplements, expands, and reinforces the ana- 
lytical methods that were, and continue to be, the main focus of this book. The analytical methods 
are the foundation of electromagnetics and provide understanding and physical interpretation of 
electromagnetic phenomena and interactions. Although numerical and computational methods 
have, especially in the last four decades, played a key role in the solution of complex elec- 
tromagnetic problems, they are highly dependent on fundamental principles. Not understanding 
the basic fundamentals of electromagnetics, represented by analytical methods, may lead to the 
lack of physical realization, interpretation and verification of simulated results. In fact, there are 
a plethora of personal and commercial codes that are now available, and they are expanding 
very rapidly. Users are now highly dependent on these codes, and we seem to lose focus on the 
interpretation and physical realization of the simulated results because, possibly, of the lack of 
understanding of fundamental principles. There are numerous books that address numerical and 
computational methods, and this author did not want to repeat what is already available in the lit- 
erature, especially with space limitations. Only the moment method (MM), in support of Integral 
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Equations (IEs), and Diffraction Theory (GTD/UTD) are included in this book. However, to aid in 
the computation, simulation and animation of results based on analytical formulations included 
in this book, even provide some of the data in graphical form, forty-eight basic MATLAB® 
computer programs have been developed and are included in the website that is part of this book. 

The first edition was based on material taught on a yearly basis and notes developed over nearly 
20 years. This second edition, based on an additional 20 years of teaching and development of 
notes and multimedia (for a total of over 40 years of teaching), refined any shortcomings of 
the first edition and added: a new chapter, two new complete sections, numerous smaller inserts, 
examples, numerous end-of-chapter problems, and Multimedia (including PPT notes, MATLAB® 
computer programs for computations, simulations, visualization, and animation). The four Fortran 
programs from the first edition were translated in MATLAB®, and numerous additional ones 
were developed only in MATLAB®. These are spread throughout Chapters 4 through 14. The 
revision of the book also took into account suggestions of nearly 20 reviewers selected by 
the publisher, some of whom are identified and acknowledged based on their approval. The 
multicolor PowerPoint (PPT) notes, over 4,200 viewgraphs, can be used as ready-made lectures 
so that instructors will not have to labor at developing their own notes. Instructors also have the 
option to add PPT viewgraphs of their own or delete any that do not fit their class objectives. 

The book can be used for at least a two-semester sequence in Electromagnetics, beyond an 
introduction to basic undergraduate EM. Although the first part of the book in intended for senior 
undergraduates and beginning graduates in electrical engineering and physics, the later chapters 
are targeted for advanced graduate students and practicing engineers and scientists. The majority 
of Chapters 1 through 10 can be covered in the first semester, and most of Chapter 11 through 15 
can be covered in the second semester. To cover all of the material in the proposed time frame 
would be, in many instances, an ambitious task. However, sufficient topics have been included 
to make the text complete and to allow instructors the flexibility to emphasize, de-emphasize, or 
omit sections and/or chapters. Some chapters can be omitted without loss of continuity. 

The discussion presumes that the student has general knowledge of vector analysis, differential 
and integral calculus, and electromagnetics either from at least an introductory undergraduate 
electrical engineering or physics course. Mathematical techniques required for understanding 
some advanced topics, mostly in the later chapters, are incorporated in the individual chapters or 
are included as appendixes. 

Like the first edition, this second edition is a thorough and detailed student-oriented book. The 
analytical detail, rigor, and thoroughness allow many of the topics to be traced to their origin, 
and they are presented in sufficient detail so that the students, and even the instructors, will 
follow the analytical developments. In addition to the coverage of traditional classical topics, 
the book includes state of the art advanced topics on DNG Metamaterials, Artificial Impedance 
Surfaces (AIS, EBG, PBG, HIS, AMC, PMC), Integral Equations (IE), Moment Method (MM), 
Geometrical and Uniform Theory of Diffraction (GTD/UTD) for PEC and impedance surfaces, 
and Green’s functions. Electromagnetic theorems, as applied to the solution of boundary-value 
problems, are also included and discussed. 

The material is presented in a methodical, sequential, and unified manner, and each chapter is 
subdivided into sections or subsections whose individual headings clearly identify the topics dis- 
cussed, examined, or illustrated. The examples and end-of-chapter problems have been designed 
to illustrate basic principles and to challenge the knowledge of the student. An exhaustive list of 
references is included at the end of each chapter to allow the interested reader to trace each topic. 
A number of appendixes of mathematical identities and special functions, some represented also 
in tabular and graphical forms, are included to aid the student in the solution of the examples 
and assigned end-of-chapter problems. A solutions manual for all end-of-chapter problems is 
available exclusively to instructors. 
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In Chapter 1, the book covers classical topics on Maxwell’s equations, constitutive param- 
eters and relations, circuit relations, boundary conditions, and power and energy relations. The 
electrical properties of matter for both direct-current and alternating-current, including an update 
on superconductivity, are covered in Chapter 2. The wave equation and its solution in rectan- 
gular, cylindrical and spherical coordinates are discussed in Chapter 3. Electromagnetic wave 
propagation and polarization is introduced in Chapter 4. Reflection and transmission at normal 
and oblique incidences are considered in Chapter 5, along with depolarization of the wave due 
to reflection and transmission and an introduction to metamaterials (especially those with neg- 
ative index of refraction, referred to as double negative, DNG). Chapter 6 covers the auxiliary 
vector potentials and their use toward the construction of solutions for radiation and scattering 
problems. The theorems of duality, uniqueness, image, reciprocity, reaction, volume and surface 
equivalences, induction, and physical and physical optics equivalents are introduced and applied 
in Chapter 7. Rectangular cross section waveguides and cavities, including dielectric slabs, arti- 
ficial impedance surfaces (AIS) [also referred to as Electromagnetic Band-Gap (EBG) structures; 
Photonic Band-Gap (PBG) structures; High Impedance Surfaces (HIS), Artificial Magnetic Con- 
ductors (AMC), Perfect Magnetic Conductors (PMC)], striplines and microstrips, are discussed in 
Chapter 8. Waveguides and cavities with circular cross section, including the fiber optics cable, 
are examined in Chapter 9, and those of spherical geometry are introduced in Chapter 10. Scat- 
tering by strips, plates, circular cylinders, wedges, and spheres is analyzed in Chapter 1 1 . Chapter 
12 covers the basics and applications of Integral Equations (IE) and Moment Method (MM). The 
techniques and applications of the Geometrical and Uniform Theory of Diffraction (GTD/UTD) 
are introduced and discussed in Chapter 13. The PEC GTD/UTD techniques of Chapter 13 
are extended in the new Chapter 14 to wedges with impedance surfaces, utilizing Maliuzhinets 
functions. The classic topic of Green’s functions is introduced and applied in Chapter 15. 

Throughout the book an e jr '“ time convention is assumed, and it is suppressed in almost all the 
chapters. The International System of Units, which is an expanded form of the rationalized MKS 
system, is used throughout the text. In some instances, the units of length are given in meters (or 
centimeters) and feet (or inches). Numbers in parentheses ( ) refer to equations, whereas those 
in brackets [ ] refer to references. For emphasis, the most important equations, once they are 
derived, are boxed. 

I would like to acknowledge the invaluable suggestions from those that contributed to the 
first edition, too numerous to mention here. Their names and contributions are stated in the first 
edition. It is a pleasure to acknowledge the invaluable suggestions and constructive criticisms of 
the reviewers of this edition who allowed their names to be identified (in alphabetical order): 
Prof. James Breakall, Penn State University; Prof. Yinchao Chen, University of South Carolina; 
Prof. Ramakrishna Janaswamy, University of Massachusetts, Amherst; Prof. Ahmed Kishk, Uni- 
versity of Mississippi; Prof. Duncan McFarlane, University of Texas at Dallas; Prof. Jeffrey Mills, 
formerly of IIT, Chicago; Prof. James Richie, Marquette University; Prof. Yahya Rahmat-Samii, 
UCLA (including additional correspondence on the topics of DNG metamaterials and Artificial 
Impedance Surfaces); Prof. Phillip E. Serafim, Northeastern University; Prof. Ahmed Sharkawy, 
University of Delaware; and Prof. James West, Oklahoma State University. There have been other 
reviewers and contributors to this edition. In addition, I would like to thank Dr. Timothy Griesser, 
Agilent Technologies, for allowing me to use material from his PhD dissertation at Arizona State 
University for the new Chapter 14; Prof. Sergey N. Makarov, Worcester Polytechnic Institute, 
for providing MATLAB® programs for computations and animations of scattering by cylinders 
and spheres; Prof. Nathan Newman, Arizona State University, for updates on superconductivity; 
Prof. Donald R. Wilton, University of Houston, for elucidations on the topic of integral equations; 
Dr. Arthur D. Yaghjian for his review and comments on the Veselago planar lens; Prof. Danilo 
Ericcolo, University of Illinois at Chicago, for bringing to my attention some updates to the first 
edition; and Dr. Lesley A. Polka, Intel, for allowing me to use figures from her MS Thesis and 
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Time-Varying and Time-Harmonic 
Electromagnetic Fields 



1.1 INTRODUCTION 

Electromagnetic field theory is a discipline concerned with the study of charges, at rest and in 
motion, that produce currents and electric-magnetic fields. It is, therefore, fundamental to the 
study of electrical engineering and physics and indispensable to the understanding, design, and 
operation of many practical systems using antennas, scattering, microwave circuits and devices, 
radio-frequency and optical communications, wireless communications, broadcasting, geosciences 
and remote sensing, radar, radio astronomy, quantum electronics, solid-state circuits and devices, 
electromechanical energy conversion, and even computers. Circuit theory, a required area in the 
study of electrical engineering, is a special case of electromagnetic theory, and it is valid when 
the physical dimensions of the circuit are small compared to the wavelength. Circuit concepts, 
which deal primarily with lumped elements, must be modified to include distributed elements and 
coupling phenomena in studies of advanced systems. For example, signal propagation, distortion, 
and coupling in microstrip lines used in the design of sophisticated systems (such as computers and 
electronic packages of integrated circuits) can be properly accounted for only by understanding 
the electromagnetic field interactions associated with them. 

The study of electromagnetics includes both theoretical and applied concepts. The theoretical 
concepts are described by a set of basic laws formulated primarily through experiments conducted 
during the nineteenth century by many scientists — Faraday, Ampere, Gauss, Lenz, Coulomb, 
Volta, and others. They were then combined into a consistent set of vector equations by Maxwell. 
These are the widely acclaimed Maxwell’s equations. The applied concepts of electromagnetics 
are formulated by applying the theoretical concepts to the design and operation of practical 
systems. 

In this chapter, we will review Maxwell’s equations (both in differential and integral forms), 
describe the relations between electromagnetic field and circuit theories, derive the boundary 
conditions associated with electric and magnetic field behavior across interfaces, relate power and 
energy concepts for electromagnetic field and circuit theories, and specialize all these equations, 
relations, conditions, concepts, and theories to the study of time-harmonic fields. 



1.2 MAXWELL’S EQUATIONS 

In general, electric and magnetic fields are vector quantities that have both magnitude and 
direction. The relations and variations of the electric and magnetic fields, charges, and cur- 
rents associated with electromagnetic waves are governed by physical laws, which are known 
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as Maxwell’s equations. These equations, as we have indicated, were arrived at mostly through 
various experiments carried out by different investigators, but they were put in their final form 
by James Clerk Maxwell, a Scottish physicist and mathematician. These equations can be written 
either in differential or in integral form. 



1.2.1 Differential Form of Maxwell’s Equations 

The differential form of Maxwell’s equations is the most widely used representation to solve 
boundary-value electromagnetic problems. It is used to describe and relate the field vectors, current 
densities, and charge densities at any point in space at any time. For these expressions to be valid, 
it is assumed that the field vectors are single-valued, bounded, continuous functions of position 
and time and exhibit continuous derivatives. Field vectors associated with electromagnetic waves 
possess these characteristics except where there exist abrupt changes in charge and current densi- 
ties. Discontinuous distributions of charges and currents usually occur at interfaces between media 
where there are discrete changes in the electrical parameters across the interface. The variations of 
the field vectors across such boundaries (interfaces) are related to the discontinuous distributions 
of charges and currents by what are usually referred to as the boundary conditions . Thus a com- 
plete description of the field vectors at any point (including discontinuities) at any time requires 
not only Maxwell’s equations in differential form but also the associated boundary conditions. 

In differential form, Maxwell’s equations can be written as 



where 
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All these field quantities — %, Dt. 92), 98, $, M, and p — are assumed to be time-varying, and 
each is a function of the space coordinates and time, that is, % — % (x,y,z; t). The definitions 
and units of the quantities are 

% — electric field intensity (volts/meter) 

Dt — magnetic field intensity (amperes/meter) 

2) = electric flux density (coulombs/square meter) 

98 = magnetic flux density (webers/square meter) 

= impressed (source) electric current density (amperes/square meter) 

$ c — conduction electric current density (amperes/square meter) 

$ d = displacement electric current density (amperes/square meter) 

Mj — impressed (source) magnetic current density (volts/square meter) 

M d — displacement magnetic current density (volts/square meter) 

^ = electric charge density (coulombs/cubic meter) 
j* = magnetic charge density (webers/cubic meter) 
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The electric displacement current density $ d — was introduced by Maxwell to complete 

Ampere’s law for statics, V x = $. For free space, $ d was viewed as a motion of bound 
charges moving in “ether,” an ideal weightless fluid pervading all space. Since ether proved to be 
undetectable and its concept was not totally reasonable with the special theory of relativity, it has 
since been disregarded. Instead, for dielectrics, part of the displacement current density has been 
viewed as a motion of bound charges creating a true current. Because of this, it is convenient to 
consider, even in free space, the entire 32)/3t term as a displacement current density. 

Because of the symmetry of Maxwell’s equations, the 328/3 t term in (1-1) has been des- 
ignated as a magnetic displacement current density. In addition, impressed (source) magnetic 
current density Mj and magnetic charge density ^ have been introduced, respectively, in 
(1-1) and (1-4) through the “generalized” current concept. Although we have been accustomed to 
viewing magnetic charges and impressed magnetic current densities as not being physically real- 
izable, they have been introduced to balance Maxwell’s equations. Equivalent magnetic charges 
and currents will be introduced in later chapters to represent physical problems. In addition, 
impressed magnetic current densities, like impressed electric current densities, can be considered 
energy sources that generate fields whose expressions can be written in terms of these current 
densities. For some electromagnetic problems, their solution can often be aided by the introduc- 
tion of “equivalent” impressed electric and magnetic current densities. The importance of both 
will become more obvious to the reader as solutions to specific electromagnetic boundary-value 
problems are considered in later chapters. However, to give the reader an early glimpse of the 
importance and interpretation of the electric and magnetic current densities, let us consider two 
familiar circuit examples. 

In Figure 1-la, an electric current source is connected in series to a resistor and a parallel- 
plate capacitor. The electric current density can be viewed as the current source that generates 
the conduction current density $ c through the resistor and the displacement current density $ d 
through the dielectric material of the capacitor. In Figure 1-1 b, a voltage source is connected to a 
wire that, in turn, is wrapped around a high-permeability magnetic core. The voltage source can 
be viewed as the impressed magnetic current density that generates the displacement magnetic 
current density through the magnetic material of the core. 

In addition to the four Maxwell’s equations, there is another equation that relates the variations 
of the current density $ ic and the charge density y . Although not an independent relation, this 
equation is referred to as the continuity equation because it relates the net flow of current out of 
a small volume (in the limit, a point) to the rate of decrease of charge. It takes the form 




The continuity equation 1-6 can be derived from Maxwell’s equations as given by (1-1) through 
(l-5c). 

1.2.2 Integral Form of Maxwell’s Equations 

The integral form of Maxwell’s equations describes the relations of the field vectors, charge 
densities, and current densities over an extended region of space. They have limited applications, 
and they are usually utilized only to solve electromagnetic boundary-value problems that possess 
complete symmetry (such as rectangular, cylindrical, spherical, etc., symmetries). However, the 
fields and their derivatives in question do not need to possess continuous distributions . 

The integral form of Maxwell’s equations can be derived from its differential form by utilizing 
the Stokes' and divergence theorems. For any arbitrary vector A, Stokes’ theorem states that the 
line integral of the vector A along a closed path C is equal to the integral of the dot product of 
the curl of the vector A with the normal to the surface S that has the contour C as its boundary . 
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Figure 1-1 Circuits with electric and magnetic current densities, (a) Electric current density, (b) Magnetic 
current density. 



In equation form, Stokes’ theorem can be written as 






( V x A) • ds 



(1-7) 



The divergence theorem states that, for any arbitrary vector A, the closed surface integral of the 
normal component of vector A over a surface S is equal to the volume integral of the divergence of 
A over the volume V enclosed by S . In mathematical form, the divergence theorem is stated as 



j (j A.- ds — JJJ VA dv (1-8) 

Taking the surface integral of both sides of (1-1), we can write 

JJ (V x *). ds = - JJ Mi • ds - JJ ^ • ds = - JJ jH, • ds - ^ JJ a,- • ds (1-9) 

Applying Stokes’ theorem, as given by (1-7), on the left side of (1-9) reduces it to 

J %-dt = - Wi-ds (l-9a) 
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which is referred to as Maxwell’s equation in integral form as derived from Faraday’s law. In the 
absence of an impressed magnetic current density, Faraday’s law states that the electromotive 
force (emf) appearing at the open-circuited terminals of a loop is equal to the time rate of decrease 
of magnetic flux linking the loop. 

Using a similar procedure, we can show that the corresponding integral form of (1-2) can be 
written as 

j> x-dl = jj s $ ic • ds + JJs a • ds = !J s $ IC ■ ds + JJ $ d ■ ds (1-10) 



which is usually referred to as Maxwell’s equation in integral form as derived from Ampere ’ s law. 
Ampere ’s law states that the line integral of the magnetic field over a closed path is equal to the 
current enclosed. 

The other two Maxwell equations in integral form can be obtained from the corresponding 
differential forms, using the following procedure. First take the volume integral of both sides of 
(1-3); that is, 
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where % is the total electric charge. Applying the divergence theorem, as given by (1-8), on the 
left side of (1-11) reduces it to 




2)- ds — 
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( 1 - 1 1 a) 



which is usually referred to as Maxwell’s electric field equation in integral form as derived from 
Gauss’s law. Gauss’s law for the electric field states that the total electric flux through a closed 
surface is equal to the total charge enclosed. 

In a similar manner, the integral form of (1-4) is given in terms of the total magnetic charge 
a m by 



ds = a„ 



( 1 - 12 ) 



which is usually referred to as Maxwell’s magnetic field equation in integral form as derived from 
Gauss’s law. Even though magnetic charge does not exist in nature, it is used as an equivalent 
to represent physical problems. The corresponding integral form of the continuity equation, as 
given by (1-6) in differential form, can be written as 




= - 



9 

dt 





(1-13) 



Maxwell’s equations in differential and integral form are summarized and listed in Table 1-1. 



1 .3 CONSTITUTIVE PARAMETERS AND RELATIONS 

Materials contain charged particles, and when these materials are subjected to electromagnetic 
fields, their charged particles interact with the electromagnetic field vectors, producing currents 
and modifying the electromagnetic wave propagation in these media compared to that in free 
space. A more complete discussion of this is in Chapter 2. To account on a macroscopic scale for 
the presence and behavior of these charged particles, without introducing them in a microscopic 
lattice structure, we give a set of three expressions relating the electromagnetic field vectors. 
These expressions are referred to as the constitutive relations, and they will be developed in 
more detail in Chapter 2. 
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TABLE 1-1 Maxwell’s equations and the continuity equation in differential and integral forms for 
time-varying fields 

Differential form Integral form 
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One of the constitutive relations relates in the time domain the electric flux density 2) to the 
electric field intensity % by 



2) = £ * % 



(1-14) 



where e is the time-varying permittivity of the medium (farads/meter) and * indicates convolution. 
For free space 



e = E 0 = 8.854xl(T 12 



1CT 9 
36 jt 



(farads/meter) 



(1-1 4a) 



and (1-14) reduces to a product. 

Another relation equates in the time domain the magnetic flux density 98 to the magnetic field 
intensity M by 

9 8 = (1-15) 



where fi is the time-varying permeability of the medium (henries/meter). For free space 

jl = // 0 = 4tt x 1 0 7 (henries/meter) (l-15a) 



and (1-15) reduces to a product. 

Finally, the conduction current density $ c is related in the time domain to the electric field 
intensity % by 

$ c =o*% (1-16) 



where a is the time-varying conductivity of the medium (siemens/meter). For free space 

rr=0 (l-16a) 

In the frequency domain or for frequency nonvarying constitutive parameters, the relations (1-14), 
(1-15) and (1-16) reduce to products. For simplicity of notation, they will be indicated everywhere 
from now on as products, and the caret Q in the time-varying constitutive parameters will be 
omitted. 

Whereas (1-14), (1-15), and (1-16) are referred to as the constitutive relations, e,jl and a are 
referred to as the constitutive parameters, which are, in general, functions of the applied field 
strength, the position within the medium, the direction of the applied field, and the frequency of 
operation. 
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The constitutive parameters are used to characterize the electrical properties of a material. 
In general, materials are characterized as dielectrics (insulators), magnetics, and conductors, 
depending on whether polarization (electric displacement current density), magnetization (mag- 
netic displacement current density), or conduction (conduction current density) is the predominant 
phenomenon. Another class of material is made up of semiconductors , which bridge the gap 
between dielectrics and conductors where neither displacement nor conduction currents are, in 
general, predominant. In addition, materials are classified as linear versus nonlinear, homoge- 
neous versus nonhomogeneous (inhomogeneous), isotropic versus nonisotropic (anisotropic), and 
dispersive versus nondispersive , according to their lattice structure and behavior. All these types 
of materials will be discussed in detail in Chapter 2. 

If all the constitutive parameters of a given medium are not functions of the applied field 
strength, the material is known as linear-, otherwise it is nonlinear. Media whose constitu- 
tive parameters are not functions of position are known as homogeneous-, otherwise they are 
referred to as nonhomogeneous (inhomogeneous) . Isotropic materials are those whose constitu- 
tive parameters are not functions of direction of the applied field; otherwise they are designated 
as nonisotropic (anisotropic). Crystals are one form of anisotropic material. Material whose con- 
stitutive parameters are functions of frequency are referred to as dispersive-, otherwise they are 
known as nondispersive. All materials used in our everyday life exhibit some degree of disper- 
sion, although the variations for some may be negligible and for others significant. More details 
concerning the development of the constitutive parameters can be found in Chapter 2. 



1.4 CIRCUIT-FIELD RELATIONS 

The differential and integral forms of Maxwell’s equations were presented, respectively, in 
Sections 1.2.1 and 1.2.2. These relations are usually referred to as field equations, since the 
quantities appearing in them are all field quantities. Maxwell's equations can also be written in 
terms of what are usually referred to as circuit quantities ; the corresponding forms are denoted 
circuit equations . The circuit equations are introduced in circuit theory texts, and they are special 
cases of the more general field equations. 



1.4.1 Kirchhoff’s Voltage Law 



According to Maxwell’s equation l-9a, the left side represents the sum voltage drops (use the 
convention where positive voltage begins at the start of the path) along a closed path C, which 
can be written as 

%-dl (volts) (1-17) 

The right side of (l-9a) must also have the same units (volts) as its left side. Thus, in the absence 
of impressed magnetic current densities (M,,- = 0), the right side of (l-9a) can be written as 
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because by definition i// m = L s i where L s is an inductance (assumed to be constant) and i is the 
associated current. Using (1-17) and ( 1 - 17a), we can write (l-9a) with M,, — 0 as 
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(1-1 7b) 



Equation l-17b states that the voltage drops along a closed path of a circuit are equal to the time 
rate of change of the magnetic flux passing through the surface enclosed by the closed path, or 
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equal to the voltage drop across an inductor L s that is used to represent the stray inductance of 
the circuit. This is the well-known Kirchhoff loop voltage law, which is used widely in circuit 
theory, and its form represents a circuit relation. Thus we can write the following field and circuit 
relations: 



(f %-dl = (( 
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Field Relation 
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dt 
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Circuit Relation 
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( 1 - 17c) 



In lumped-element circuit analysis, where usually the wavelength is very large (or the dimen- 
sions of the total circuit are small compared to the wavelength) and the stray inductance of the 
circuit is very small, the right side of (1-1 7b) is very small and it is usually set equal to zero. In 
these cases, (l-17b) states that the voltage drops (or rises) along a closed path are equal to zero, 
and it represents a widely used relation to electrical engineers and many physicists. 

To demonstrate Kirchhoff s loop voltage law, let us consider the circuit of Figure 1-2 where 
a voltage source and three ideal lumped elements (a resistance R, an inductor L, and a capacitor 
C) are connected in series to form a closed loop. According to (1-1 7b) 



di 

— V s + Vr + Vl + Vc — —L s — — —V s i 

dt 



(1-18) 



where L s , shown dashed in Figure 1-2, represents the total stray inductance associated with the 
current and the magnetic flux generated by the loop that connects the ideal lumped elements (we 
assume that the wire resistance is negligible). If the stray inductance L s of the circuit and the 
time rate of change of the current is small (the case for low-frequency applications), the right 
side of (1-18) is small and can be set equal to zero. 



1.4.2 Kirchhoff’s Current Law 



The left side of the integral form of the continuity equation, as given by (1-13), can be written 
in circuit form as 




(1-19) 



+ V R 




Figure 1-2 RLC series network. 
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where ^ i represents the sum of the currents passing through closed surface S. Using (1-19) 
reduces (1-13) to 
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(C s v) = 




(l-19a) 



since by definition % = C s v where C s is a capacitance (assumed to be constant) and v is the 
associated voltage. 

Equation l-19a states that the sum of the currents crossing a surface that encloses a circuit 
is equal to the time rate of change of the total electric charge enclosed by the surface, or equal 
to the current flowing through a capacitor C s that is used to represent the stray capacitance of 
the circuit. This is the well-known Kirchhojf node current law, which is widely used in circuit 
theory, and its form represents a circuit relation. Thus, we can write the following field and circuit 
relations: 



Field Relation 



Circuit Relation 



JJ' , d fff , d% d% dv 

L t “- ds = -V,jjj v ^ iv = * Z' = -nr = “ C -T7 



dt 



dt 



(1-1 9b) 



In lumped-element circuit analysis, where the stray capacitance associated with the circuit is very 
small, the right side of (l-19a) is very small and it is usually set equal to zero. In these cases, 
(l-19a) states that the currents exiting (or entering) a surface enclosing a circuit are equal to zero. 
This represents a widely used relation to electrical engineers and many physicists. 

To demonstrate Kirchhoff’ s node current law, let us consider the circuit of Figure 1-3 where 
a current source and three ideal lumped elements (a resistance R, an inductor L, and a capacitor 
C) are connected in parallel to form a node. According to (1-1 9a) 



dv 

+ ir + if, + ic = C s — — = — i s c (1-20) 

dt 

where C s , shown dashed in Figure 1-3, represents the total stray capacitance associated with the 
circuit of Figure 1-3. If the stray capacitance C s of the circuit and the time rate of change of 
the total charge % are small (the case for low-frequency applications), the right side of (1-20) is 
small and can be set equal to zero. The current i s c associated with the stray capacitance C s also 



S 




Figure 1-3 RLC parallel network. 
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includes the displacement ( leakage ) current crossing the closed surface S of Figure 1-3 outside of 
the wires. 



1.4.3 Element Laws 

In addition to Kirchhoff’s loop voltage and node current laws as given, respectively, by (l-17b) 
and (l-19a), there are a number of current element laws that are widely used in circuit theory. 
One of the most popular is Ohm’s law for a resistor (or a conductance G), which states that the 
voltage drop Vr across a resistor R is equal to the product of the resistor R and the current iR 
flowing through it {vr — RiR or i R = vr/R — Gvr). Ohm’s law of circuit theory is a special case 
of the constitutive relition given by (1-16). Thus 



Field Relation Circuit Relation 

$ c = a% 4» i R = -v R = Gv r 
K 



( 1 - 21 ) 



Another element law is associated with an inductor L and states that the voltage drop across an 
inductor is equal to the product ofL and the time rate of change of the current through the inductor 
(vi — LdiL/dt). Before proceeding to relate the inductor’s voltage drop to the corresponding held 
relation, let us hrst define inductance. To do this we state that the magnetic flux f m is equal to the 
product of the inductance L and the corresponding current i . That is — Li . The corresponding 
held equation of this relation is (1-15). Thus 



Field Relation Circuit Relation 
S® = yu3€ f m = Li L 



( 1 - 22 ) 

Using (l-5c) and (1-15), we can write for a homogeneous and non-time-varying medium that 

(l-22a) 



a m 9 , ™ 

M.i = = — (ji at) = U 

dt dt dt 



where M d is dehned as the magnetic displacement current density [analogous to the electric 
displacement current density $ d — 39)/3f = d(e%)/dt — ed%/dt], With the aid of the right side 
of (l-9a) and the circuit relation of (1-22), we can write 



3 

dt 



ff. 



. , d'I'm 9 ,, . . J 9 «'l 

• ds — = — (Lit) — L — = vr 

dt dt 



i s 9r 

Using (l-22a) and (l-22b), we can write the following relations: 

Circuit Relation 



Field Relation 

.. dm 

'™'d — M M < ^ 
ot 



, di L 
v L — L — 
dt 



(l-22b) 



(l-22c) 



Using a similar procedure for a capacitor C, we can write the held and circuit relations 
analogous to (1-22) and (l-22c): 



Field Relation 

3 = e% 



= e at 



Circuit Relation 
4^ % — C v e 

3 vc 

* ic = c nr 



(1-23) 

(1-24) 



A summary of the held theory relations and their corresponding circuit concepts are listed in 
Table 1-2. 



TABLE 1-2 Relations between electromagnetic Held and circuit theories 
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1.5 BOUNDARY CONDITIONS 

As previously stated, the differential form of Maxwell’s equations are used to solve for the field 
vectors provided the field quantities are single-valued, bounded, and possess (along with their 
derivatives) continuous distributions. Along boundaries where the media involved exhibit discon- 
tinuities in electrical properties (or there exist sources along these boundaries), the field vectors 
are also discontinuous and their behavior across the boundaries is governed by the boundary 
conditions . 

Maxwell’s equations in differential form represent derivatives, with respect to the space coor- 
dinates, of the field vectors. At points of discontinuity in the field vectors, the derivatives of the 
field vectors have no meaning and cannot be properly used to define the behavior of the field 
vectors across these boundaries. Instead, the behavior of the field vectors across discontinuous 
boundaries must be handled by examining the field vectors themselves and not their derivatives. 
The dependence of the field vectors on the electrical properties of the media along boundaries of 
discontinuity is manifested in our everyday life. It has been observed that cell phone, radio, or 
television reception deteriorates or even ceases as we move from outside to inside an enclosure 
(such as a tunnel or a well-shielded building). The reduction or loss of the signal is governed 
not only by the attenuation as the signal/wave travels through the medium, but also by its behav- 
ior across the discontinuous interfaces. Maxwell’s equations in integral form provide the most 
convenient formulation for derivation of the boundary conditions. 

1.5.1 Finite Conductivity Media 

Initially, let us consider an interface between two media, as shown in Figure 1-4 a, along which 
there are no charges or sources. These conditions are satisfied provided that neither of the two 
media is a perfect conductor or that actual sources are not placed there. Media 1 and 2 are 
characterized, respectively, by the constitutive parameters e\,p.\,o\ and Si, P- 2,02 ■ 

At a given point along the interface, let us choose a rectangular box whose boundary is denoted 
by Co and its area by So- The x, y, z coordinate system is chosen to represent the local geometry 
of the rectangle. Applying Maxwell’s equation l-9a, with Mi — 0, on the rectangle along Co and 
on So, we have 

<j> %-dl = ~ jfjf 28-r/s (1-25) 

As the height Ay of the rectangle becomes progressively shorter, the area So also becomes 
vanishingly smaller so that the contributions of the surface integral in (1-25) are negligible. In 
addition, the contributions of the line integral in (1-25) along Ay are also minimal, so that in the 
limit (Ay — > 0), (1-25) reduces to 



-a x Ax — ? 2 -a,Ar = 0 

%i t — %n = 0 =>■ %\ t — %i t (1-26) 



or 



n x (?2 — ^ 1 ) = 0 



01,02 are finite 



(l-26a) 



In (1-26), %\ t and %n represent, respectively, the tangential components of the electric field in 
media 1 and 2 along the interface. Both (1-26) and (l-26a) state that the tangential components 
of the electric field across an interface between two media, with no impressed magnetic current 
densities along the boundary of the interface, are continuous . 
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(b) 

Figure 1-4 Geometry for boundary conditions of tangential and normal components, (a) Tangential. 
(b) Normal. 



Using a similar procedure on the same rectangle but for (1-10), assuming = 0, we can write 
that 

2£tr - W 2t = 0 =► = W 2t (1-27) 



or 



n x (3€i — 3Ci) = 0 



<ti,CT 2 are hnite 



(l-27a) 



which state that the tangential components of the magnetic field across an interface between two 
media, neither of which is a perfect conductor, are continuous . This relation also holds if either 
or both media possess hnite conductivity. Equations l-26a and l-27a must be modified if either 
of the two media is a perfect conductor or if there are impressed (source) current densities along 
the interface. This will be done in the pages that follow. 

In addition to the boundary conditions on the tangential components of the electric and mag- 
netic fields across an interface, their normal components are also related. To derive these relations, 
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let us consider the geometry of Figure 1-4 b where a cylindrical pillbox is chosen at a given point 
along the interface. If there are no charges along the interface, which is the case when there are 
no sources or either of the two media is not a perfect conductor, (1-1 la) reduces to 



d s = 0 






(1-28) 



As the height Ay of the pillbox becomes progressively shorter, the total circumferential area A\ 
also becomes vanishingly smaller, so that the contributions to the surface integral of (1-28) by 
A i are negligible. Thus (1-28) can be written, in the limit (Ay — ► 0), as 



®2-a y A 0 - 2)i - a y A 0 = 0 

2>2„ - % n = 0 =► 3 2n = 2) ln (1-29) 



or 



n • (9> 2 — 9>i) = 0 



(Ti, cr? are finite 



(l-29a) 



In (1-29), 2)i„ and 2 ) 2 „ represent, respectively, the normal components of the electric flux density 
in media 1 and 2 along the interface. Both (1-29) and (l-29a) state that the normal components 
of the electric flux density across an interface between two media, both of which are imperfect 
electric conductors and where there are no sources, are continuous. This relation also holds if 
either or both media possess finite conductivity. Equation l-29a must be modified if either of the 
media is a perfect conductor or if there are sources along the interface. This will be done in the 
pages that follow. 

In terms of the electric field intensities, (1-29) and (l-29a) can be written as 



£2&2n — £\&ln = 4 - %2 n — — %\n => %\ n — — %2n 
£2 £l 



n- (£ 2^2 — £ 1 ^ 1 ) — 0 



o ' \ , o 2 are finite 



(1-30) 

(l-30a) 



which state that the normal components of the electric field intensity across an interface are 
discontinuous . 

Using a similar procedure on the same pillbox, but for (1-12) with no charges along the 
interface, we can write that 



®2n - n — 0 =4> 2S 2 « = 9Sl„ 



(1-31) 



n-(S(2 — 95 1) = 0 



(l-31a) 



which state that the normal components of the magnetic flux density, across an interface between 
two media where there are no sources, are continuous . In terms of the magnetic field intensities, 
(1-31) and ( 1 -3 la) can be written as 



n — n 



L 2 n 



ft, 

l±2 



^1 n 



M 2„ 

v: 



^2 n 



(1-32) 



n • (p .2 Wi — Mi i ) = 0 



(l-32a) 



which state that the normal components of the magnetic field intensity across an interface are 
discontinuous . 
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1.5.2 Infinite Conductivity Media 

If actual electric sources and charges exist along the interface between the two media, or if either 
of the two media forming the interface displayed in Figure 1-4 is a perfect electric conductor 
(PEC), the boundary conditions on the tangential components of the magnetic held [stated by 
(l-27a)] and on the normal components of the electric flux density or normal components of 
the electric held intensity [stated by (l-29a) or (l-30a)] must be modihed to include the sources 
and charges or the induced linear electric current density ($ s ) and surface electric charge density 
(jr ). Similar modihcations must be made to (l-26a), (l-31a), and (l-32a) if magnetic sources 
and charges exist along the interface between the two media, or if either of the two media is a 
perfect magnetic conductor (PMC). 

To derive the appropriate boundary conditions for such cases, let us refer hrst to Figure 1- 
4 a and assume that on a very thin layer along the interface there exists an electric surface 
charge density y (C/m 2 ) and linear electric current density $ s (A/m). Applying (1-10) along the 
rectangle of Figure l-4n, we can write that 

(f X-dt = ff $ ic .ds + — ff 2b • ds (1-33) 

Jco JJso JJ So 

In the limit as the height of the rectangle is shrinking, the left side of (1-33) reduces to 

lim (f> 9C- dl = (3€, - 2C).a T Ax (l-33a) 

Jcq 

Since the electric current density $ ic is confined on a very thin layer along the interface, the first 
term on the right side of (1-33) can be written as 

lim // $ ic -d s 

JJs 0 

= lim [$ ic • a z Ax Ay] = lim [($ ic Ay)-a z Ax] = $ s a z Ax (l-33b) 

Ay->0 Ay-*0 

Since So becomes vanishingly smaller as Ay — >■ 0, the last term on the right side of (1-33) reduces 
to 

lim — ff 2 i-ds — lim — ff 2b-a^ds—0 (l-33c) 

Ay-»0 dt JJ Sq Ay~+0 dt JJ Sq 

Substituting (l-33a) through (l-33c) into (1-33), we can write it as 



(9Ki - 9€ 2 )-a JC Ax = $ s -a- A.r 



or 





(3Kt - ^2.)-a v - $ s -a z = 0 


(l-33d) 


Since 


a A - = a v x a z 


(1-34) 


( 1 -33d) can be written as 


(3€i - 3C 2 ) • (a v x aj - $ s -a z = 0 


(1-35) 


Using the vector identity 


AB xC= CA x B 


(1-36) 
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on the first term in (1-35), we can then write it as 



a z • [(*, - X 2 ) x a v J - $, - a- = 0 (1-37) 

or 

{[a y x (3€ 2 — )] — $ s }-a z — 0 (l-37a) 

Equation l-37a is satished provided 

fi, x(* 2 -*i)-* s =0 (1-38) 



or 

a y x (Sk -*,) = #, 



(l-38a) 



Similar results are obtained if the rectangles chosen are positioned in other planes. Therefore, 
we can write an expression on the boundary conditions of the tangential components of the 
magnetic field, using the geometry of Figure 1-4 a, as 



A x(X 2 -*!) = *, 



(1-39) 



Equation 1-39 states that the tangential components of the magnetic field across an interface, along 
which there exists a surface electric current density $ s (A/m), are discontinuous by an amount 
equal to the electric current density . 

If either of the two media is a perfect electric conductor (PEC), (1-39) must be reduced to 
account for the presence of the conductor. Let us assume that medium 1 in Figure 1-4 a possesses 
an infinite conductivity (eri = oo). With such conductivity % \ — 0, and (l-26a) reduces to 



nx?i=04 ^2 1 = 0 



Then (1-1) can be written as 



(1-40) 



Vx«i=0 = -=4><a|=0=4>3€i=0 (1-41) 

3 1 

provided in is finite. 

In a perfect electric conductor, its free electric charges are confined to a very thin layer on the 
surface of the conductor, forming a surface charge density ^ (with units of coulombs/square 
meter). This charge density does not include bound (polarization) charges (which contribute 
to the polarization surface charge density) that are usually found inside and on the surface of 
dielectric media and form atomic dipoles having equal and opposite charges separated by an 
assumed infinitesimal distance. Here, instead, the surface charge density p represents actual 
electric charges separated by finite dimensions from equal quantities of opposite charge. 

When the conducting surface is subjected to an applied electromagnetic field, the electric 
surface charges are subjected to electric field Lorentz forces. These charges are set in motion and 
thus create a surface electric current density $ s with units of amperes per meter. The surface 
current density $ s also resides in a vanishingly thin layer on the surface of the conductor so that 
in the limit, as Ay — > 0 in Figure 1-4 a, the volume electric current density $ (A/m 2 ) reduces to 



lim (^ Ay) = $ s 

Ay-»0 

Then the boundary condition of (1-39) reduces, using (1-41) and (1-42), to 

nxl 2 =^4 3^2 1 = $ s 



(1-42) 



(1-43) 
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which states that the tangential components of the magnetic field intensity are discontinuous next 
to a perfect electric conductor by an amount equal to the induced linear electric current density . 

The boundary conditions on the normal components of the electric field intensity, and the 
electric flux density on an interface along which a surface charge density y resides on a very 
thin layer, can be derived by applying the integrals of (1-1 la) on a cylindrical pillbox as shown 
in Figure 1-4 b. Then we can write (1-1 la) as 

lim 0 S i-ds — lim /// a dv (1-44) 

JJao.Ai Aj^O JJJ v 7ev 

Since the cylindrical surface A\ of the pillbox diminishes as Ay -> 0, its contributions to the 
surface integral vanish. Thus we can write (1-44) as 

(S 2 - ®i) -n A 0 = A lim o [( ?gu Ay)T 0 ] = <f g A 0 (1-45) 



which reduces to 

n-(2) 2 - ®i) = <j es =>• n ~ 2>i„ = y es 



( 1 -45 a) 



Equation l-45a states that the normal components of the electric flux density on an interface, 
along which a surface charge density resides, are discontinuous by an amount equal to the surface 
charge density. 

In terms of the normal components of the electric field intensity, (l-45a) can be written as 



n- (£ 2^2 — £l^l) = </ es 



(1-46) 



which also indicates that the normal components of the electric field are discontinuous across a 
boundary along which a surface charge density resides. 

If either of the media is a perfect electric conductor (PEC) (assuming that medium 1 possesses 
infinite conductivity o\ = oo), (l-45a) and (1-46) reduce, respectively, to 



A ®2 = ?es =► = ( /es 



n-^2 = y e J e 2 => %2n = </ e j£ 2 



(l-47a) 

(l-47b) 



Both (l-47a) and (l-47b) state that the normal components of the electric flux density, and corre- 
sponding electric field intensity, are discontinuous next to a perfect electric conductor. 



1.5.3 Sources Along Boundaries 

If electric and magnetic sources (charges and current densities) are present along the interface 
between the two media with neither one being a perfect conductor, the boundary conditions on 
the tangential and normal components of the fields can be written, in general form, as 



-n x (% 2 - «i) = M s 
n x (3 <e 2 - 9ti) = $ s 
M®2 - ®i) = f es 
n-(98-i — ®i) = a 



(l-48a) 

(l-48b) 

(l-48c) 

(l-48d) 



ms 
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TABLE 1-3 Boundary conditions on instantaneous electromagnetic fields 





General 


Finite 

conductivity media, 
no sources or 
charges 
< 71,02 # OO 

h = 0; 9es = 0 

M s = 0; a =0 


Medium 1 of 
infinite electric 
conductivity 
<*1 = *! = 0) 

< 7 l = OO; 02 7^ OO 

M, = 0 ; a =0 


Medium 1 of 
infinite 
magnetic 
conductivity 
(81 = = 0) 

*,=o;?„ = o 


Tangential 
electric field 
intensity 


— n x (8 2 — 8i) = M-s 


n x ('82 — 81) =0 


n x 8 2 = 0 


— n x 8 2 = it s 


Tangential 
magnetic field 
intensity 


n x(* 2 -*,) = & 


n x (3^2 — = 0 


A x ve 2 = $ s 


n x 3^2 = 0 


Normal electric 
flux density 


n- ( 9 ) 2 - 3 i) =^ ej 


n • ( 9 ) 2 — 2>!) = 0 


fl S 2 =? es 


n • S& 2 = 0 


Normal magnetic 
flux density 


n.(SL-S&,) = ?m( 


n- ( 3 2 - ®i) = 0 


n ■ 98 2 = 0 


n • ®2 = (/ 

f ms 



where and (</ m , ) are the magnetic and electric linear (per meter) current and surface 

(per square meter) charge densities, respectively. The derivation of (l-48a) and (l-48d) proceeds 
along the same lines, respectively, as the derivation of (l-48b) and (l-48c) in Section 1.5.2, but 
begins with (l-9a) and (1-12). 

A summary of the boundary conditions on all the field components is found in Table 1-3, 
which also includes the boundary conditions assuming that medium 1 is a perfect magnetic 
conductor (PMC). In general, a magnetic conductor is defined as a material inside of which both 
time-varying electric and magnetic fields vanish when it is subjected to an electromagnetic field. 
The tangential components of the magnetic field also vanish next to its surface. In addition, the 
magnetic charge moves to the surface of the material and creates a magnetic current density that 
resides on a very thin layer at the surface. Although such materials do not physically exist, they 
are often used in electromagnetics to develop electrical equivalents that yield the same answers 
as the actual physical problems. PMCs can be synthesized approximately over a limited frequency 
range (band-gap); see Section 8.8 . 



1 .6 POWER AND ENERGY 

In a wireless communication system, electromagnetic fields are used to transport information over 
long distances. To accomplish this, energy must be associated with electromagnetic fields. This 
transport of energy is accomplished even in the absence of any intervening medium. 

To derive the equations that indicate that energy (and forms of it) is associated with electro- 
magnetic waves, let us consider a region V characterized by s, n, o and enclosed by the surface S, 
as shown in Figure 1-5. Within that region there exist electric and magnetic sources represented, 
respectively, by the electric and magnetic current densities and M, . The fields generated by 
and Mi that exist within S are represented by 8. Ht. These fields obey Maxwell’s equations, 
and we can write using (1-1) and (1-2) that 



Vx? = — Mj = — Mj — /x — = —Mi — M d 

3 1 3 1 



(l-49a) 
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Figure 1-5 Electric and magnetic fields within S generated by and M, . 



V x X = $, : + $ c + — = $, + + e— = + $ c + $ d (l-49b) 

dt dt 

Scalar multiplying (l-49a) by 3€ and (l-49b) by IS, we can write that 

9€-(V x %) = + M d ) (l-50a) 

*.(Vx3£) = «.(l,+l c +W ( 1 -50b) 

Subtracting (l-50b) from (l-50a) reduces to 

* • ( V x *) - * • (V x *) = -31 C-(M, +M d ) -%■($, +$ c +$ d ) (1-51) 

Using the vector identity 

V • (A x B) = B • ( V x A) — A • (V x B ) (1-52) 

on the left side of (1-51), we can write that 

V.<*xX) = -X.(M i +JL d )-%.& i +$ c +} d ) (1-53) 



or 



V • (« X 3K) + +JH d )+9- (J, + $ c + $ d ) = 0 



Integrating (1-53) over the volume V leads to 

30 dv — — HI ^ ‘(^i + ^d) + ^ • (<?, + <?f + $d)]d\) 
Applying the divergence theorem (1-8) on the left side of (1-54) reduces it to 

%xW).ds = - [[[ [3 te.(M,+M d ) + %.($! +$ c + $ d )]dv 
S JJJ V 

or 




(l-53a) 



(1-54) 



(1-55) 



(l-55a) 



Equations l-53a and l-55a can be interpreted, respectively, as the differential and integral forms 
of the conservation of energy . To accomplish this, let us consider each of the terms included in 
(l-55a). 
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The integrand, in the first term of (l-55a), has the form 



$P = %x% e (1-56) 

where £P is known as the Poynting vector. It has the units of power density (watts/square meter), 
since % has units of volts/meter and ’Hi has units of ampere/meter, so that the units of are 
volts ■ ampere/meter 2 = watts/meter 2 . Thus the first term of (l-55a), written as 



< 3>e = jf (* x X) ■ ds = JJ Sf-ds (1-57) 

represents the total power 2P e exiting the volume V bounded by the surface S . 

The other terms in (l-55a), which represent the integrand of the volume integral, can be written 
as 



where 



90S 93€ 1 9X 2 9/1 9 

3 t-Jid = 3€ = /z3C = -n = — -/zX 2 = —w n 

dt dt 2^ dt dt \2^ ) dt 

/, d =%.$ c = %.(a%) = a% 2 

„ „ 9 Si) 9* 1 d% 2 9 /I ,\ 9 

d dt dt 2 dt dt 1 2 / dt 



(l-58a) 

(l-58b) 

(l-58c) 

(l-58d) 



m m — -/zX 2 — magnetic energy density (J/m 3 ) (l-58e) 

w e — -e/' 2 = electric energy density(J/m 3 ) (l-58f) 

A = — supplied power density(W/m 3 ) (l-58g) 

A — fT ’/ 2 = dissipated power density(W/m 3 ) (l-58h) 



Integrating each of the terms in (l-58a) through (l-58d), we can write the corresponding forms 
as 



S ?s = - JJJ (X-Mi+t-Mdv = JJJ /> s dv (l-59a) 

III <*■■*'> <•’ = l III = b III dv = (1 - 59b) 

S ? d = JJJ (%-$ c )dv = JJJ o% 2 dv = JJJ /> d dv (l-59c) 

fflw-bffl k^-bllfy-^-b^ (1 - 59d) 



